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Abstract 

By applying a magnetic field whose Zeeman energy exceeds the Kondo energy by 
an order of magnitude the ground state of the Friedel- Anderson impurity is a magnetic 
state. In recent years the author introduced the Friedel Artificially Inserted Resonance 
(FAIR) method to investigate impurity properties. Within this FAIR approach the 
magnetic ground state is derived. Its full excitation spectrum and the composition of 
the excitations is calculated and numerically evaluated. From the excitation spectrum 
the electron density of states is calculated. Majority and minority d-resonances are 
obtained. The width of the resonances is about twice as wide as the mean field theory 
predicts. This broadening is due to the fact that any change of the occupation of the 
d-state in one spin band changes the eigenstates in the opposite spin band and causes 
transitions in both spin bands. This broadening reduces the height of the resonance 
curve and therefore the density of states by a factor of two. This yields an intuitive 
understanding for a previous result of the FAIR approach that the critical value of the 
Coulomb interaction for the formation of a magnetic moment is twice as large as the 
mean field theory predicts. 

PACS: 75.20.Hr, 71.23.An, 71.27.+a 
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1 Introduction 



The properties of magnetic impurities in a metal is one of the most intensively studied prob- 
lems in solid state physics. The work of Friedel [1] and Anderson [2j laid the foundation to 
understand why some transition-metal impurities form a local magnetic moment while oth- 
ers don't. Kondo [3] showed that multiple scattering of conduction electrons by a magnetic 
impurity yields a divergent contribution to the resistance in perturbation theory. Yoshida 
[4J introduced the concept that the (spin 1/2) magnetic impurity forms a singlet state with 
the conduction electrons and is non-magnetic at zero temperature. These new insights stim- 
ulated a large body of theoretical and experimental work (see for example [1], [5], [6], [7], 

0, m)- 

The majority of experimental and theoretical work has focussed on the singlet Kondo 
ground state. However, the "magnetic state" of the impurity is of equal or even greater 
importance because magnetic impurities are always present, including in micro-chips and 
nanostructures, and influence the thermodynamic and transport properties of the hosts. 
Since many experiments and almost all technical applications are not performed at low 
temperatures the magnetic impurities are generally far above their Kondo temperature 
and show their full magnetic behavior. The theoretical investigation of the magnetic state 
has been explored in much less detail than the Kondo ground state for spin 1/2 impurities. 

In many cases the Kondo temperature is very low, in the range of liquid helium tem- 
perature. In this case the impurity is in the magnetic state at relatively low temperature. 
(The word impurity is in this paper reserved to impurities which possess - at sufficiently 
high temperature - a magnetic moment). When the temperature is several times the Kondo 
temperature one is sufficiently above to destroy the Kondo ground state. On the other 
hand one may expect that the properties of the magnetic state are not yet influenced by 
the thermal excitations due to the finite temperature. Therefore a number of theoretical 
investigations treat the magnetic state at zero temperature, i.e. as a magnetic ground state. 
This approach is probably justified but it leaves the work always vulnerable to the criticism 
that there is no magnetic moment at zero temperature. 

Therefore in this paper I prefer to use the effect of a magnetic field on the Kondo state. 
A magnetic field which is an order of magnitude larger than fcsTx/yUB (yUB=Bohr magneton) 
destroys the Kondo singlet state as well and yields the magnetic state. Its side effects are 
that it changes the energy of the d-states by iyU^i? and shifts the conduction bands by 
±^bB. The latter yields the Pauli susceptibility but has otherwise only a negligible effect 
on the interaction between the impurity and the conduction electrons because the Fermi 
level for spin-up and down electrons readjusts to the same height (as before). 

Friedel [1] and Anderson [2] derived a criterion for the instability of the paramagnetic 
state, i.e. the formation of a magnetic moment: Take the density of states N(i{sf) of the 
d-resonance at the Fermi energy (in the paramagnetic state) and multiply it by the Coulomb 
repulsion energy U . If the product NdU > 1 then a magnetic moment is formed. Within 
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mean field theory the d-density of states is given by a Lorentz function 
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where Ed^a is an effective energy of the d-electrons in the spin-up or down state, Ed^a = 
Ed + U (ud-a) while {rid -a) is the average occupation of the d-electron with the opposite 
spin) and the resonance width F^/ is given in mean field theory by 

Here Vsd is the s-d-hopping matrix element between a conduction electron and the d-state 
at the impurity and Ng is the density of states of the conduction electrons. In the mean field 
theory an occupied electron state can only make transitions into cj^^-states. (Throughout 
this paper I express electron states by their creation operators). 

It is well known that the mean field theory has a number of shortcomings. During the last 
few years the group of the author has developed a new approach to the impurity problem, 
in particular the Friedel- Anderson and the Kondo impurity. In this FAIR method a Friedel 
state is Artificially built from each conduction band and Inserted as a Resonance state into 
the conduction or s-band of spin-up and spin-down electrons. In the appendix a short review 
of the FAIR solution for the Friedel impurity is sketched. 

The FAIR solution for the magnetic state yields a considerably lower energy for the 
"magnetic ground state" and requires a much larger critical Coulomb interaction to form a 
magnetic state. This is of some practical importance since the mean field approximation is 
used in a number of numerical spin-density functional theory calculations for the magnetic 
moment of impurities in an (s,p) metal host [2], [15], [1^], [I7|, [IB]. 

In addition to the size of the magnetic moment one would like to know the density of 
states in the magnetic state. The answer of the mean field theory has been discussed above. 
But there have been a number of suggestions that the d-resonance is broader than the mean 
field suggests (see for example Logan p^j). The mean field theory decouples the spin-up 
d-electron from the spin-down d-electron, but in reality the d-electrons are coupled through 
the Coulomb energy. A transition in the (i| electron state changes the energy and the state 

of the (i| electron as well. Therefore it has been suggested in the past that the d-resonances 
in the Friedel-Anderson impurity are larger than the mean field theory predicts. A wider 
d-resonance in the Friedel-Anderson impurity together with the condition NdU > 1 would 
require a larger Coulomb energy to form a magnetic moment. In this connection the previous 
result of the author that the FAIR solution requires a (two times) larger Coulomb energy to 
form a magnetic moment would find a simple physical interpretation. 

It is the goal of this paper to calculate the density of states of the "magnetic ground state" 
in the FAIR solution and compare it with the mean field density of states. In section II the 
theoretical background of the magnetic state of the Friedel-Anderson impurity is sketched. 
In section III electrons and holes are introduced into the magnetic ground state. Their 
interactions and excitation energies are derived. In section IV the results of the numerical 
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calculations are presented. Finally in section V and VI the results are discussed together 
with the conclusion. In the appendix A the basic idea of the FAIR method is sketched. 



2 Theoretical Background 

The simplified Hamiltonian for a magnetic impurity is generally described by the Friedel- 
Anderson (FA) Hamiltonian 

{N-1 N-1 ^ 

J2 ^AaCva + EddU^ + Vf[dlc,„ + \ + UUd^Udl (1) 

i/=0 v=0 ) 

Here the operators cl^ represent s-electrons, i.e. the conduction band. 
2.1 The FAIR method 

In the Friedel-Anderson Hamiltonian in equ. ([1]) the d-state for each spin interacts with 
every electron in the conduction band. Imagine how much easier the task would be if the 
d-electron would interact only with a single electron state (in each spin band). All other 
conduction band states would represent just a background or quasi- vacuum. This is the 
FAIR approach. 

During the last few years the author introduced such a solution to the Friedel-Anderson 
impurity problem in which only four electron states, the spin-up and spin-down d-states (i| 
and (i| and two FAIR states, and interact through the Coulomb and s-d-hopping 

potential. These states aj^ and are composed of the spin-up and spin-down conduction 
band states. They are the Friedel Artificially Inserted Resonance states or FAIR states. The 
interaction of the remaining conduction electron states with the d-states is insignificant; they 
just yield a background. This yields very good ground-state properties. The FAIR states 
are composed of the corresponding conduction bands 

EN-l y 7 s:^N-lnv 

The remaining (A^ — 1) states in each spin band are constructed orthogonal to the corre- 
sponding FAIR state, orthonormal to each other and sub-diagonal with respect to the band 
energy Hamiltonian 

N-1 
i/=0 

This yields new bases for the conduction bands and with 1 < i < (A^ — 1). 

These new bases are uniquely determined by the two FAIR states. 
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Within the new bases the FA-Hamiltonian ([T]) can be expressed as 



with 



where 



N-l 



i=l 

N-l N-l 



(2) 
(3) 



i=l 



i=l 



and the spin-down Hamiltonians are obtained by replacing t by | and the aj-states by 
6|-states. 

2.1.1 Nest states 

The Hamiltonian H'q is diagonal in the band states a\-^ and h\ ^ for < i < — 1. The 
only interaction takes place between the states aj ^ rf-f, 6o,j, and d^. I call these states the 
nest states. The ground state of the Hamiltonian Hq is straight forward. It consists of the 
coupled state between the nest states and a partially occupied spin-up and down band. I 
occupy each spin component with N/2 electrons, putting n = N/2 — 1 electrons into each 
conduction band states and one spin- up and one spin-down electron into the nest. This 
yields the magnetic ground state as described in equ. (jl]). 



^Ms= Aa,bal^bl^ + Aa,dal^dl + Ad,bd\bl^ + Ad,dd\d\ \0a,^0b,i) (4) 

where |0a,^0b,j) = 11^=1 '^]t nj=i \^o) represents a kind of quasi-vacuum {n = N/2). 

The calculation of the coefficients Aa,b, ■■ yields a secular Hamiltonian Hy^ which 1 call 
the nest-Hamiltonian and which has the form 



Tjnst 
^1/1 





ysd 


ysd 

a 





ysd 


E't^ + Ed 





ysd 


ysd 
a 





Ed + E^^ 


ysd 





ysd 
a 


ysd 


2Ed + U 



Here the abbreviations are used: V^'^ = V^^"-'*'^ and Vf^'^ = Vq^^^'^ (see equ. ([2])). The superscript 
nst stands for nest and the subscript 1/1 gives the number of nest electrons in the spin-up 
and spin-down state. (The energy of the occupied band states is not included. It yields 
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the same contribution to each component). Hy^ has four eigenvalues and eigenstates. The 
lowest eigenvalue yields the ground state. The components of the ground state, nest plus 
the band states, are shown in Fig.l. 

The first order correction to the energy, i.e. the expectation value of H[ is zero. But in 
addition the second order perturbation of H[ is extremely small. This is demonstrated in 
appendix B. 

It may appear remarkable that the neglect of the interactions between the d-electron and 
all the band states |oj || and jfej || yields a realistic ground state. But it is not unheard off 
that one can obtain an excellent ground state while neglecting a major part of the interaction 
in the system. The BCS theory is a good example because it only includes the electron- 
phonon interaction between Cooper pairs of time-reversed electrons. The interaction between 
all the other electrons is neglected although their number is much larger. 

One major part of the numerical calculation is, of course, the optimization of the two 
FAIR states and 6q so that the expectation value of the energy Eqq = (^ms l-f^ol ^^ms) 
of the Hamiltonian Hq has a minimum. The optimization procedure is described at length 
in previous papers [20], [21], [22] and is taken for granted in this paper and will not be 
described here. (The FAIR states are rotated in Hilbert space). Since we don't count the 
FAIR states any more as band states the number of band states is reduced by one and their 
energy is slightly shifted (by less than the original energy spacing). The band states enter 
in the energy Eqq only through the kinetic (band) energy of the occupied band states. In a 

way they just prepare the nest for the states aj |, dj, 6o,j,' 
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Fig.l: The composition of the magnetic state \1/ms in the nest is shown. 
It consists of four Slater states. Each Slater state has a half-full spin-up and 
down band, two FAIR states (circles in within the bands) and two d-states 
(circles on the left and right of the band). Full black circles represent 
occupied states and light grey represent empty states. The band at the right 
with the half-filled circles symbolizes the magnetic solution with four 
Slater states. 



In the numerical calculation we will present the results for two examples with the pa- 
rameters U = 1.0, Ed = —0.5 and |^ = 0.05 and |^ = 0.025. The smaller value of 
the s-d-matrix element permits a better fit of the resulting resonance curve with a Lorentz 
curve since the effect of the finite width of the band is smaller. 
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2.1.2 Self-consistent perturbation 

In the construction of the magnetic ground state \1/m5 the Hamiltonian H[ has been com- 
pletely neglected. Below we will derive the excitation energies by introducing an additional 
electron (hole) into an empty (occupied) states. For this calculation it is important to know 
whether the empty state is really empty or whether transitions from the ground state into 
the state due to H[ have partially occupied this state. (This problem is well known from 
the calculation of the electro n-phonon mass enhancement. In the calculation of the electron- 
phonon self-energy one injects an electron into an "empty state" k above the Fermi energy. 
The transitions of this electron via the electron-phonon interaction into other empty states 
k' contribute to the self-energy S. However, the state k was not really empty because tran- 
sitions from the ground state into k already created a finite occupation of k. One has to 
correct the self-energy due to these processes). 

In the appendix I show that transitions from the ground state into empty band states 
(due to H[) are practically zero. The interference between transitions from the d-state and 
from the corresponding FAIR state almost perfectly cancel each other. The total weight in all 
perturbation states is only of the order of 10~^ and can be completely neglected. Therefore 
the band states are either completely empty or fully occupied. 



3 Calculation of Excitations 
3.1 Injection of an electron 

In Fig. 2 a spin-up electron is injected into one of the empty states aj^ of the |a||-band. 

This yields the Slater states (a). The Slater state (/5) is obtained by injecting an electron 
into the nest, either into the state aL or 
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Fig.2: An electron has been injected into the spin-up band and the spin up nest. 
This induces electron or hole transitions in the spin-down band. The resulting 
Slater states are shown as (7) and (6). (Each band with half circles consists 
of two Slater states.) 
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The injection into or (i| yields 

al^'^Ms = (^Ad,bald\bl^ + Ad4ald]d\^ \Oa,'iOb,i) for 

Both final states yield a double occupancy of the spin up nest states. Furthermore these 
states are not eigenstates of the nest. With respect to the (basis) states ajc/jfe,!,^ |0a,|06,|) and 

|0a,T0b,j.) the nest Hamiltonian takes the form 



sd 
b 



By diagonalization one obtains the eigenstates \E'2/\ and \l/2/\ with 

vl/(;J = at^4(<)6t^ + 5(f4)|0,,T0,,) (7) 

The states (a) and (/?) in Fig.2 are the initial states which one obtains through injection 
of a spin-up electron into the ground state. Due to the perturbation Hamiltonian H[ these 
states interact with each other and interacts with the states (7) and {5). 

In table I the possible states which can be obtained through the injection of a spin-up 
electron plus linear coupling through H[ are collected. These states are (a) the two-electron 
nest ground state plus one electron, {(3) a nest with two spin- up and one spin-down electron, 
(7) a full spin-up and empty spin-down nest plus one spin down electron and (5) a full spin- 
up and full spin-down nest and one spin down hole. In table I these states, their number 
and their energies are listed. 





number 


energy 






N/2 






^2/1 


2, a = 1,2 


Ed + -E'o"^ + -^2/1 ~ -^00 






N/2 


Ed + E\^^ + - Eoo 






N/2-1 


2Ed + U + E^"^ + E^o^ 





Table I: This table describes the states of Fig.2 which are generated 
by the injection of one spin-up electron into the magnetic ground state 
and transition from the resulting states through H[. The state '^ms is 
given by equ. (j4]) and \1/2"| is given by ([7]). The 
energy is measured from the ground-state energy Eqq. 

Fig.2 and table I show all the spin-up electron excitations which interact linearly in H[. 
This is a total of (|A^ + l) states. It is straight forward to construct the secular matrix (i.e. 
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the excitation Hamiltonian H^'^*) between the excitations in table I. One may put the two 
nest states ^^2/1 ^2/1 positions one and two, followed by the (|A^ — l) additional 

single particle excitations. We denote these (|A^ + l) states as ip^. The diagonal of the 
Hamiltonian is given by the energies in table I. The off-diagonal elements of H^'^^ are the 
matrix elements of H[ between the states {a, f3, 7, 6) in Fig.2. The single particle excitations 
interact only with the first two nest states through H[ but not among each other. (In 
appendix B the corresponding matrix elements are shown in table III for similar transitions 
from the ground state.) 

This Hamiltonian H^'^* is diagonalizcd and yields a set of (|A^ + l) new eigenstates 
with cigenenergies EJ^^^. The components of the eigenstates ipfj, in terms of Lp,, are given as 
columns ('(/'^) . The z/th row of the matrix [ip^) yields the amplitude of our z/th original state 
(fi, in terms of the new eigenstates ip/^. If this uth original state is, for example, aj^^'MS then 
it can be expressed in the new eigenstates as 



Its density of states is then 



Nuie) = EM" S {e - Ef) 



Since electron injection creates only the states (a) and one obtains the full (spin-up) 
excitation spectrum by summing over these (^^N + 2) states. The weight of states at^^^^s 
is one, however, the weight of aj^^jv^s is only |^d,6|^ + since 



Adfiol-fdlbli + Ad,dal^d\d\ \ lOa^^Ob,;) 



A similar result is found for the weight of dj^MS which is |>lo,b|^ + |^a,(i|^- 

Since Oq^^'ms and (ij^MS are not eigenstatct of Hq they represent a combination of the 
two eigenstates ^2/1 • Fig-3 is sketched what happens when an electron is injected into 
either the state or dj. The electron injection yields a superposition of the two eigenstates 

"^yl- From these states the electron can make a transition into any of the (7) states via 
H[. The two transition amplitudes interfere in this transition. This interference has to 
be included in the calculation of the spectral weight density of a^^^^MS and (i|\t'MS (which 

requires just the scalar products between al^'^MS and ^'2"^ (or rfj^'Ms and ^2/1)- '^^^^ 
discussed in more detail in appendix C. 
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Fig.3: An electron has been injected into the aj^ state (or cij state, 
dashed arrow). The resulting state is a superposition of two nest states. 
From these nest states the electron makes (as one possibility) a transition 
into the state a^^ where the two amplitudes interfere. The energies of the 
two nest states are shown. The different thickness of the arrows shows 
different probabilities for the two paths. 



3.2 Injection of a hole 

For the full spectrum of excitations one has to include the injection of holes into the occupied 
states. This is shown in Fig. 4. The hole can be injected into the occupied states a^^ yielding 
aj^'^MS or into the nest. In the latter case the spin-up part of the nest is emptied. This 
yields for the secular matrix of the nest in analogy to equ. 

Hsn = ( ^j. \ ) (8) 
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Fig.4: An hole has been injected into the spin-up band and the spin up nest. 
This induces electron or hole transitions in the spin-down band. The resulting 
Slater states are shown as (7) and (5). 
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state 


number 


energy 




N/2 - 1 






2, a = 1,2 






N/2 - 1 






N/2 





Table II: This table describes the states of Fig. 4 which are generated 
by the injection of one spin-up hole into the magnetic ground state 
and transition from the resulting states through H[. The energy is 
measured from the ground-state energy £'oo- 



The construction of the excitation or secular Hamiltonian if^'^* is in complete analogy to 
the electron injection. This time the number of excitations is |iV. The spectrum is obtained 
in the same way as before. 

4 Numerical Results 

For the numerical calculation a conduction band with a finite number of states is used. We 
follow here Wilson [23] by using an s-electron band with constant density of states and the 
Fermi level in the center, which we divide into energy cells (Liy. In each energy cell (which 
may contain k-states c|^) we rearrange the states (by an orthogonal transformation) so 
that one state cj, = Zu^^'^Yle '-k accumulates all the interaction with the d-states while the 
other {Z^ — 1) states have zero interaction with the d-states. Wilson normalized the energy 
in terms of the Fermi energy so that his band extended from —1 to 1. Wilson's logarithmic 
scale for the energy cells is not opportune for the present investigation because it is not 
fine enough at the energy of the d-resonance. Therefore I use a linear sub-division of the 
energy band (—1 : 1). For the majority of calculations the energy band is sub-divided into 

= 40, 80 and 160 energy cells. The state cj, represents all the s-electron states in the cell 

= (— 1 + : —1 + + I) j^) and possesses the average energy = —1 + (^u + ^) 
(corresponding to —.975, —.925, .. -|- .975 for N = 40). For energy cells with constant 
width of 2/iV the s-d- matrix elements VJ'^ is given by V^^/Vn so that \ VJ'^\'^ = iKdl^- 

In the following I show the results for the Friedel- Anderson Hamiltonian with the d-level 
energy Ed = —0.5, the Coulomb energy U = 1.0 and an s-d-hopping matrix element of 
iKdl^ = 0.025. The magnetic moment of this impurity is /i = 0.998/iB (for A^ = 80). The 
calculations are performed with A^ = 40, 80 or 160 energy levels of constant spacing. A 
second set of results is derived for the parameters = —0.5, U = 1 and \V^j'^ = 0.05. 

To assure that the magnetic state "^ms is indeed the ground state (i.e. to prevent 
the formation of a singlet state) a magnetic field B is applied yielding a magnetic energy 
Eb = f^sB. This energy is chosen so that Eb > lOksTx- I estimate the Kondo energy from 
the difference between the energies of the singlet and the triplet state. This energy difference 
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is about 8 x 10^^ for |^ = 0.05 and of the order of 10"^ for \V^°f = 0.025. It turned 
out that the required magnetic energy is in both cases so small that it yields no noticeable 
changes. This is partly due to the fact that the absolute smallest band energies which are 
±1/N act as a finite temperature as Wilson pointed out p3]. For = 40 this corresponds 
to a temperature of £i?/40 which is a very large temperature compared with most Kondo 
temperatures. But magnetic field is of academic importance to assure the magnetic state 
"^Ms is the appropriate ground state. 

The use of equidistant energy levels is important to identify the resonance state within 
the electron bands. But it has the drawback that it does not describe well the behavior of 
the wave function at low energies. At low energies the logarithmic energy scale which Wilson 
introduced would be more appropriate. But the evaluation of the density of states is much 
more difficult for a non-linear energy scale. 

For each spin band one obtains a spectrum with a total weight of (A^ + 1), corresponding 
to A^ s-electron states and one d-electron state. However, the number of energy levels is 
(3A^ + 1). (This is the number of eigenstates of the excitation Hamiltonians for electrons 
and holes together). This means that the weight at the individual energies is at least for 
2A^ energies much less than one. In Fig.5 the spectral weight at different energies is shown 
in the energy range from —1 to +1 for the minority band. For negative energies the weight 
is either very close to one or very small. Here one can calculate the density of states by 
the separation of the levels with weight close to one (by dividing the weight by the level 
separation). The evaluation is more complicated for positive energies. However, here the 
sum of neighboring energy levels is close to one. Then one can calculated the "center of 
weight" for two neighboring levels which have a total weight close to one and then proceed 
as before. It turns out that the best approach is to start from the lower and upper ends of 
the band in the evaluation. 
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Fig.5: The spectral weight for the different excitation energies. 

In Fig. 6a the density of states of the excitation spectrum for the majority spin is shown. 
The full circles are obtained with N — AO states and the stars use N — SO equally spaced 
Wilson states. The full curve represents a Lorentz curve with the parameters 

, , 1 0.08 

Nd £ = 5 

7r(£- (-.53))' + 0.082 

The resonance energy is is Er = —0.53 and the resonance (half) width is = 0.08. The 
corresponding mean field resonance width is F^/ = 0.039. 
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energy 



Fig. 6a: The density of states for the majority spins. The full circles 
are calculated with iV = 40 states and the stars use N = 80 equally 
spaced Wilson states. The full curve represents a Lorentz curve 
with the resonance energy Er = —0.53 and the width = 0.08. 
This width is twice the mean-field value of F^/ = 0.039. 

In Fig. 6b the density of states of the minority spin is drawn. Again the full curve 
represents a Lorentz curve with the resonance at Er = 0.52 and a resonance half- width of 
F^ = 0.8. 
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Fig. 6b: The density of states for the minority spins. The full circles 
are calculated with iV = 40 states and the stars use = 80 equally 
spaced Wilson states. The full curve represents a Lorentz curve 
with the resonance energy Er = 0.52 and the width = 0.08. 
Again this width is twice the mean-field value of F^/ = 0.039. 

For comparison Fig. 7 shows the result of a similar calculation and evaluation for a simple 
Friedel resonance, where the Coulomb energy U is set equal to zero. The evaluation yields 
a Lorentz curve with = —0.51 and F^ = 0.044. The s-d-matrix element is still given by 
WPf = 0.025. 
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energy 

Fig. 7: The density of states for an impurity with U — 0. The full circles 
and the stars are calculated with = 40 and N = 80 states. 
The full curve represents a Lorentz curve with the resonance energy 
Er = -0.51 and the width = 0.044. 

Alternatively we tried to obtain the density of states by broadening the 5-shaped energy 
spectrum in Fig. 5 with a Gaussian curve ^l/27rexp [{e - E^f /2a'^] . This method worked 
quite well. The optimal density curve was obtained when \/2a was equal to the level distance 
2/N . The width of the two resonances was essentially the same as in Fig.6a,b and Fig.8a,b. 
Only the heights were slightly reduced. However, I prefer to use the other evaluation method 
in this paper so that there is no doubt that the broadening of the resonance width is a real 
physical effect and not due to an artificial broadening with a Gaussian curve. 

In a second series the majority and minority density of states have been calculated for 
the parameters E^ = —0.5, U = 1, and \Vg^f = 0.05. This impurity possesses a magnetic 
moment oi /i — 0.997 fj,B- In this case the mean-field theory yields a resonance width of 
rm/=0.079. The best fit to the numerical results yields Tfair — 0.17. It is again about 
twice the value of the mean field. (The resonance curves are no longer perfectly symmetrical 
because of the finite width of the conduction band). 

The corresponding Friedel density of states (which is not shown here) has a resonance 
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with of 0.08 which is quite close to the theoretical value of 0.079. 



2.0n 




energy 

Fig. 8a: The density of states for the majority spins for \V^X = 0.05 



2.0n 




energy 

Fig. 8b: The density of states for the minority spins for \V^/ = 0.05 
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5 Discussion 



In the mean field approximation the magnetic state has two d-resonances at the energies 
Ed^fj = Ed + Uud^^a- Since in the symmetric case one has n^^o- = (1 =f /z) /2 one finds 
Ed^fj = Ed + U {1^ jj) /2 = =i=f/ lijl = ±fiEd since Ed + U/2 = 0. So generally the resonances 
are closer to Fermi energy than ±Ed. In our case for the parameters Ed = —0.5, U = 1 and 
11^,°^! = 0.025 the magnetic moment is within 1% /i » 1 (in units of /i^) and one expects the 
resonance almost at ±0.5. The resonance width in mean field is given by T^f = vr | V"^^| g (ep) 
where g (ep) is the density of states of the s-electrons at the Fermi level. This yields for the 
above parameters F^/ = vr * 0.025 * .5 = 0.039. (The matrix element V^^ and the density of 
states are normalized to the atomic volume as the sample volume). This is the same width 
that one expects for a Friedel resonance with IV;^^! = 0.025. Indeed in Fig.7 the Friedel 
resonance has a F- value of F^r = 0.044. This agrees within 10% or 0.005 with the numerical 
result. 

The important result is that the resonance width of the minority and majority spins is 
larger than the Friedel resonance width by a factor of two. Therefore the obtained resonance 
width is also twice the mean field resonance width. This suggests that the any calculation 
which uses mean field yields an incorrect density of states. It will be interesting to check 
how the spin-density functional theory is affected by this result because the latter uses the 
mean field approximation. 

Since the Coulomb interaction broadens the d-resonance by a factor of two it also reduces 
the height of the resonance by the same factor of two. Therefore it is very plausible that the 
mean field theory overestimates the tendency to form a magnetic moment. If the criterion 
for the formation of a magnetic moment, UNd > 1, is accepted then one expects that the 
critical Coulomb energy for the formation of a magnetic moment is increased by a factor of 
two. This was the previous result by the author [2T]. 

In the density of states of the majority and minority spins in Fig.Gab and Fig.8a,b one 
observes a scattering and a small maximum at zero energy. This is probably due to the fact 
that I used a constant cell width for the Wilson states. This means that I average over all 
states within an energy cell of the width 6E = 2/N. This is definitely a poor approximation 
for the two energy cells Cn/2-i and (which touch the Fermi level). Wilson avoids this 
problem by using a logarithmic energy scale. However, the present method to evaluate the 
density of states does not work for an energy-dependent cell width. Details of this question 
will have to be clarified in the future. 

6 Conclusion 

In this paper the density of states of the Friedel-Anderson impurity is calculated in the 
magnetic ground state. The magnetic ground state is enforced by the application of a 
magnetic field whose Zeeman energy is an order of magnitude larger than the Kondo energy. 
(For the parameters chosen in the numerical calculation the effect of the magnetic field is so 
small that it can be neglected). The FAIR ground state is the eigenstate of a Hamiltonian 
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Hq. To construct this Hamiltonian two artificial Priedel resonance states aj^ and bl ^ are 
reverse engineered out of the original spin-up and spin-down conduction bands. There is 
a numerical procedure to optimize these two FAIR states. When this is done one has an 
extremely simple Hamiltonian Hq. In Hq only the two s-states aj ^ and 6q ^ interact with the 

d-state. These states a^p b^ p (i| and rfj are called the nest states. Hq is diagonal in the 
(modified) conduction band electrons |a]|| and j|- 

A perturbation Hamiltonian H[ = HpA — Hq which is the difference between the original 
FA-Hamiltonian and Hq, has zero energy expectation value in the FAIR ground state. In 
addition it is shown in second order self-consistent perturbation theory (see appendix C) 
that the total occupation of all perturbation states is only of the order of 10"^, i.e. the FAIR 
ground state has still an amplitude of 0.9998. This is important when an electron or hole is 
injected into the ground state. 

The excitation spectrum is obtained by injecting an electron or a hole into the ground 
state. The resulting excited states interact via the perturbation Hamiltonian H[ and yield a 
spectrum of energy resonances. It turns out that the injection of an electron into the spin- 
up conduction band creates also transition between the spin-down conduction band and the 
nest states. 

The resulting density of states possesses the shape of a resonance curve. However, the 
resonance width is about twice the value of the mean-field theory. As a consequence the 
height of the resonance density of states is reduced by a factor of two. Since the formation 
of a magnetic moment depends on the product of the Coulomb interaction and the density 
of d-states, one would expect that the mean field overestimates the tendency towards a 
magnetic moment. Indeed I observed in the first paper about the magnetic ground state 
that the formation of a magnetic moment requires about twice the Coulomb energy that 
the mean field theory predicts. This consorts well with the present finding of the reduced 
resonance density of states. 

The next step in the future investigation is the calculation of the density of states of 
the Kondo resonance within the FAIR model. For this calculation one has to use a Wilson 
spectrum with a logarithmic energy scale. 
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A The FAIR approach 



A.l The Friedel impurity 

The basic idea of the FAIR method can be best explained for a Priedel resonance with the 
Hamiltonian 



N-l 



N-1 



1^=0 



u=0 



This is done in the following steps: 



1. from the free (or s-) electron basis an artificial Priedel resonance (FAIR) state is 
constructed with aj = '^^QqcI together with a full orthonormal basis ^(^^ so that 
the free electron-electron Hamiltonian Hfe = ^vC^Cy takes the form 



JV-l 



N-l 



ia)fr 



i=l 



i=l 



The requirement that the matrix elements between different and a|, {i, A; ^ 0) vanish 
has the consequence that a given FAIR state aj determines uniquely the full basis |<^l|- 
In the new basis the total Friedel Hamiltonian takes the form 

Hp — Hq + H[ 

where 

7V-1 



(9) 



i=l 



The perturbation Hamiltonian has the form 

'iV-l N-l 



(10) 



i=l 



i=l 



Here the new matrix elements are given as 
2. A trial state ^i;' is defined as 



Aoal + Add''] \0a) 



where \0a) — YYj=i ^4 1^°) , n — N/2. The right side is abbreviated as 
where, for example, for a = o? one has = d^ \0a) 



(11) 
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3. The energy expectation value Eqq = {^f \Hf \ ^f) of the Hq with respect to the trial 
state \E'i? is calculated. The contribution of H[ is not included. This results in a 2x2 
secular matrix 



whose lowest eigenvalue yields Eqo and the corresponding eigenvector yields the coef- 
ficients AQ,Ad. 

4. The FAIR state is rotated (by variation) in the A^-dimensional Hilbert space until the 
lowest eigenvalue of the secular matrix reaches a minimum. 

It has been shown by the author |2l], [25] that this procedure results in the exact n- 
particle ground state of a Friedel Hamiltonian (given by ffTTl) ). The matrix elements of the 
perturbation Hamiltonian H[ between this ground state and any excited state vanish. 

It is interesting to look at the result in some more detail. 

• The states a] {i ^ 0) enter the secular matrix only through the total energy of the 
occupied states and contribute only the the background energy. 

• The coefficients and the relative weight of the states a\ and S are only determined by 
the energies of the FAIR and the d-state e'^\ Ed and their coupling V^^"'*'*'^. 

In a way one can say that the states a\ prepare just the background - a kind of nest - 
for Oq and d) . The secular matrix represents an effective Hamiltonian for these two states in 
the nest. In the following I will call the secular matrix without the kinetic energy ^21=1-^1 
the nest Hamiltonian. 



The state Oq represents an artificially inserted Friedel resonance state. Therefore I call aj 
a "Friedel Artificially Inserted Resonance" state or FAIR-state. The use of the FAIR-states 
is at the heart of my approach to the FA- and Kondo impurity problem. Therefore I call 
this approach the FAIR method. 

A. 2 From mean field to the FAIR magnetic state 

The Hamiltonian of the Friedel- Anderson impurity is given in equ. ([1]). One obtains the 
mean-field Hamiltonian from equ.(??) by replacing n^indj. =>nd^ ('^dj.) + i^di) — {nd-[) {ndi)- 
After adjusting (udi) and (udi) self-consistently one obtains two Friedel resonance Hamilto- 
nians with a spin-dependent energy of the da-state: Ed^a = Ed + U {ud-a}- 






Cua + Edadlda + ^ Vsdil^)[dlCua + cl^d^] 



N 
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The mean-field wave function is a product of two Friedel ground states for spin up and down 
Now we express each Friedel ground state "^p^ by the FAIR solution, for example 



n-1 



i=l 



For the two Friedel states in the mean-field wave function I use the form of equ. flTTl) 
and obtain for the mean-field solution 



mf 



n-1 



1=1 

Aa,b(4)iboi + Aa^ddoid^^ + Ad,bd\bl^ + Ad,dd\d^f^ |0a,|0b,j) 



n-1 



i=l 



(12) 



where and two (different) bases of the A^-dimensional Hilbert space. This 

solution can be rewritten as equation (jlj). 

In the mean-field solution "^mf the coefficients ^a,/? are restricted by two conditions 
-|. + i = 1 (A^ i -|- = 1). Therefore this state does not describe well the correlation 

effects. 

In contrast the state (jl]) opens a wide playing field for improving the solution: (i) The 



FAIR states and 6q can be individually optimized, each one defining a whole basis 

and |^]|- This yields a much better treatment of the correlation effects. The resulting 

state is denoted as the (potentially) magnetic state "^ms- The magnetic state "^ms has the 
same structure as the mean field solution \E'm/; the only difference is that its components are 
optimized for the Friedel- Anderson Hamiltonian. 



B Self-consistent Perturbation 



In the construction of the magnetic ground state "^ms only the Hamiltonian ifg (equ. 
|2]) has been used. The expectation value of the "perturbation" Hamiltonian H[ is zero, 
{■^MS \H[\ ^Ms) = where H[ = + and 



N-1 



N-1 



(13) 



But H[ yields transitions from the ground state "^ms into excited states. From equ. (fT3!) 
one recognizes that only permits transitions between the nest ^dj, aj^^ and a band state 

a^l but no transition among band states. 

If one considers only final states through this transitions which are linear in H[ then 
Fig. 9 shows the possible final states. For the spin-up band this are the transitions from the 
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nest state into an electron excitation Oj, leaving the spin- up part of the nest empty or a 
transition from an occupied state by creating a hole aj and filling both state Qq-^ and rfj of 
the spin- up part of the nest. The corresponding states are generated for the spin-down band. 




Fig. 9: Final states which can be obtained by transitions from the 
ground state "^ms through the perturbation Hamiltonian H[. 

As an example we consider the transitions in the spin-up band. Applying to the 
ground state yields 



N-1 



y{a)sd^ 



d,d 



+ 
+ 



+ vl^'>''Aa,d) a, AAA 



} |0a,T06a) 



The top two lines represent an electron excitation, shown as the left state in Fig. 9 and 
the bottom two lines a hole excitation shown as the second state in Fig.9. The final states 
and the corresponding matrix elements are collected in table III. 









excitation 


^MS 




[vj^^^'A^^. + vl^^-'A,,,) 


electron, j > n 




«j,T4 |0a,T06,i) 


{V}''^^'A,,, + V}^>'A,,,) 


electron, j > n 


^MS 




(-v;.('^^^^A,,, + \//'^)^%,) 


hole, < j < 77, 






{-V}''^^'A,,, + vj''^''A^,,) 


hole, < j < n 



Table III: For the spin-up band the final states \I'/, and the matrix elements 
of H[ between the ground state \I/ms and the final states are listed. 
The final states still have to be expanded into new eigenstates of the nest 
plus band. For the spin-down band one obtains equivalent transitions. 
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There are two important aspects of this result, (i) there are always two transitions into 



each final state, for example a, 



t 



aji and dl 



These two transitions interfere and 



cancel each other almost completely (as will be shown below), (ii) The first two final states 
in tabic III (as well as the third and forth final state) are not eigenstatcs of H^. One has to 
expand these final states in terms of the eigenstates of the nest, for example 



■J,T*2/1 



(2) 



This is in complete analogy to the calculation of the excitations in section III. 

In Fig. 10 the logarithm of the effective matrix element is plotted for the first transition 
in table III. This represents an electron excitation which is restricted to positive energies. 
As one recognizes the effective matrix elements are strongly reduced in the energy range 
in which transitions are possible. The value of Vg// lies in the range between 10^'^ and 
10^^ while the original matrix elements for Vgd sue 0.025. This applies for all eight possible 
excitations. In the energy range where an excitation is permitted the matrix elements are 
strongly reduced. Among the eight possible transitions there is only one transition whose 
matrix elements exceed 10"^. This is the second one in table III where the matrix element 
reaches values of 2 x 10~^, still much smaller than the original matrix elements. 



-1.5- 



-2.0- 



Ui -2.5- 
O 



-3.0- 



-3.5- 



Org181_8a 



-1.0 



-0.5 



magnetic state 

|V,r = 0.025 



0.5 



d 

U = 1 
N =40 



0.0 0.5 
energy 



1.0 



Fig. 10: The logarithm of the effective matrix elements for a transition 
from the ground state "^ms into the excited state Oj-^bl^ \Oa,iOb,i)- This 
transition is only possible for positive energies, and there the 
interference between the the d-state rfj and the FAIR state 
cancels the transition almost completely. 
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The strong reduction of the matrix elements is due to the introduction of the FAIR 
states which compensate the transitions involving the d-states. It is also the reason why the 
magnetic ground state is so well represented by "^ms- 

With the eigenenergies of the states in Fig. 9 and the matrix elements for the transition 
from the ground state ^'ms into these states one can now perform a self-consistent pertur- 
bation calculation. The number of excited states is 2{N — 1) for each spin-band. One can 
build the the full secular matrix for these states which consists of 4 (A'" — 1) -|- 1 states, where 
the additional one is the ground state. The result for our standard example {Ed = —0.5, 
U = 1, 1^ = 0.025, = 40) is rather dramatic. The weight of all 4 (A^ — 1) excitations 
together is only about 10""^. The amplitude of ^'m5 after diagonalization is 0.99985. This 
demonstrates the FAIR solution is an excellent approximation to the ground state. 



C Green's functions 

With the secular Hamiltonian H^'^* in section 3.1 one can construct the Green's functions 
(GF)of the excitations (p^ 

or 

G = {e + ts - 

The resulting diagonal elements of G are the Green's function of the excitations, for example 
Gil {^) is the Green's function of ^'2/1- Since the G^^^ (e) are a function of the energy the 
above relation is not very practical for a numerical nor analytical calculation of G^^^ (e). For 
a numerical evaluation it is much easier to calculate the eigenvalues and eigenvectors 
ipfj, of H^'* where ■0/* — J^u'^^Vi^- Then the Green's function 0^^^ of the excitation can be 
expressed as 



Ef,"* + is 



Each energy eigenvalue contributes to the spectrum of the state Lp,^. 

For the particle density of states we need the Green's function (or the spectrum) of the 
Y states aj^^^MS and oJ-^^ms and d|*M5- The A^/2 electron Green's functions Ga^^a^ (^) 
follow directly from the above calculation. The Gao,ao {^) and Gd,d (s) for spin up have to be 
derived from the Green's functions of "^^-^ and ^2/1- 

In Fig. 3 a spin-up electron is injected into the nest. This electron can be injected into 
the or the (i| state. These processes yield complementary amplitudes of 

Oqi^ms = (^Ad,baod\bl^ + Ad,dald\dl^ |0a,T0fe,i) for 4 
4*MS = - (Aa,bal^d\bl^ + Aa,dal^d\d\) \0a,^0b,i) for d\ 



25 



Each of the two resulting states has to be expanded in the two eigenstates "^^^ and \E'2/i 
of the nest. For example, the amplitudes for aJ|\E'A/5 are obtained through the scalar product 
between a^^^^MS and Then they can make a transition into any of the other states with 
a single excitation, for example the state o^j^'^ms- Again the resulting amplitudes interfere. 
So the particle components of the GFs Gaca,, and Gd,d (e) are composed of the GFs of 
^^/J with a weight which is the square of the scalar product between the amplitude in equ. 

([5]) and \1'2/|- These GFs Cacao (^) and Gd4 i^) have also a hole component. 
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